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Our objective is to compute the moments of the deep-inelastic structure functions of the nucleon on the lattice. 
A major source of uncertainty is the renormalization of the lattice operators that enter the calculation. In this 
talk we compare the renormalization constants of the most relevant twist-two bilinear quark operators which 
we have computed non-perturbatively and perturbatively to one loop order. Furthermore, we discuss the use of 
tadpole improved perturbation theory. 



1. Introduction 

In [1] we have initiated a lattice calculation of 
the moments of the deep-inelastic structure func- 
tions of the nucleon, both for unpolarized and 
polarized beams and targets. We will not repeat 
the results nor the details of the calculation here, 
but we refer the reader to this work. 

The program amounts to computing forward 
nucleon matrix elements of composite quark and 
gluon operators of definite twist and spin in a 
systematic fashion. Among the leading operators 
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Here {•••} means symmetrization of the indices. 
The argument /j, indicates that the operators are 
renormalized at the scale pt. 

The calculation is done in two steps. First one 
constructs bare lattice operators 0(a), which be- 
long to an irreducible representation of the hy- 
percubic group H (4) and in the limit of zero lat- 
tice spacing coincide with the classical euclidean 
continuum operators, and computes their matrix 
elements between nucleon states. The lattice op- 
erators are obtained from (2) and (4) by replacing 
the covariant derivative by the lattice covariant 
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derivative (modulo factors of i) 

These operators are in general divergent as the 
lattice spacing goes to zero. In the second step the 
bare lattice operators are renormalized by defin- 
ing 



( q (p)\0(j,)\q(p)) = ( q (p)\0(a)\ q (p))\l^ 2 
and 

{q(p)\0(n)\q(p)) = Z (n 2 ){q(p)\0(a)\q(p)), 



(6) 
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where \q(p)} is a quark state of momentum p. In 
the continuum limit this prescription amounts to 
the momentum subtraction scheme. In the fol- 
lowing we shall define the continuum quark fields 
to be \/2k times the lattice quark fields. 

In [1] we have seen that the nucleon matrix el- 
ements of the bare lattice operators can be com- 
puted relatively precisely. The renormalization 
constants Zo, on the other hand, were calcu- 
lated in lattice perturbation theory to one loop 
order. It has been argued that the perturba- 
tive series converges slowly due to the contri- 
bution of tadpole diagrams [2]. In view of this 
we have computed the renormalization constants 
non-perturbatively [3] following a recent sugges- 
tion of [4]. In this talk we shall present some 
first results of our calculation and compare non- 
perturbative and perturbative values. 

2. Non-Perturbative Renormalization 

We shall restrict ourselves to the operators 
listed in Table 1. We have nothing to add to 
the perturbative calculation. A most complete 
list of perturbative renormalization constants is 
given in [1]. For earlier results and recent inde- 
pendent work see also [5]. 

The idea of the non-perturbative calculation is 
to compute the matrix element (q(p)\0(a)\q(p)} 
in (7) directly on the lattice by numerical simu- 
lations. The only unpleasant feature in this cal- 
culation is that one has to fix the gauge. We 
have chosen the Landau gauge. This gauge suf- 
fers from the Gribov ambiguity. We sum over all 
Gribov copies. 
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Table 1 

We consider the lattice bare Green function 
Go (p) = J2 e- l(px - py) Wx)OHy))- (8) 

The operators can be written O = ipVip. We 
define the free vertex operator y tree by replac- 
ing the lattice covariant derivative by the ordi- 
nary derivative, which amounts to setting 11^ = 1 
in (5). The renormalization constant is then ob- 
tained by 



2 _ TiTS(p)V tTee (p)S(p) 



Tr TGo(p) 



Z^ 2 ) (9) 



for a suitably chosen projection matrix T. Here 
S(p) is the quark propagator and \/Z^ is the 
renormalization constant of the fermion field. 

The calculations are done on the same lattices 
as in [1]: 16 3 • 32 at (3 = 6.0 and for quenched Wil- 
son fermions with k = 0.155,0.153 and 0.1515. 
The results, using 125 configurations, are shown 
in Fig. 1 for k = 0.153. We have labelled the Z's 
by the operator matrix element they renormal- 
ize, and lvc stands for local vector current. If not 
stated otherwise, we set a = 1. 

The momenta p were chosen such that the in- 
dividual components, p^, differed by at most one 
unit of 27r/16, so that even for our largest mo- 
menta none of the p^ exceeded 7r/2. This should 
keep O(a) effects small. On the other hand, for 
the method to work, the loop momenta squared 
should be much larger than the square of a typ- 
ical hadron mass. Taking the glueball mass as a 
scale, this means that at [3 = 6.0 we can only ex- 
pect to find agreement with perturbation theory 
for 2 <C p 2 ■ The operators which involve higher 
powers of covariant derivatives are furthermore 
expected to be particularly sensitive to finite vol- 
ume effects in the region of lower momenta. 
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Figure 1. The renormalization constants Z-y yc , Za, Q , Zy 2 , Zy 3 , Za, 2 and Zy 4 as a function of fi 2 for k = 
0.1530. The solid lines are the perturbative results. The dotted and dashed lines are predictions of 
tadpole improved perturbation theory. 
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3. Comparison with Lattice Perturbation 
Theory 

Let us now compare our results with the pre- 
dictions of perturbation theory. The solid lines 
in Fig. 1 are the results of standard one-loop lat- 
tice perturbation theory [1]. There is considerable 
freedom in improving lattice perturbation theory. 
The dashed lines correspond to tadpole improve- 
ment [2], where we have taken k c = k c uo with 
k c = 0.15693 and uo = 0.891, while the dot- 
ted lines represent the choice k c = 1/8. In both 
cases the expansion parameter was taken to be 
a = 0.198 [2]. The latter option (i.e. the dotted 
line) does not fit the data any better than stan- 
dard perturbation theory. From Zy 2 onwards it 
performs even worse, so that we shall not pursue 
this choice of parameters any further. In the case 
of Zy 3 it has been found that the operator given 
in Table 1 mixes with the operator in another 
representation. The effect of mixing turned out 
to be very small though [1], so that we may safely 
neglect it here. We should also mention that for 
our choice of projection matrix T the definition 
(9) does not precisely correspond to the renor- 
malization prescription employed in perturbation 
theory. The two definitions can be converted into 
each other. The difference is insignificant. 

The figures are organized according to the 
power of covariant derivatives in the operators. 
In one-loop perturbation theory Z\ yc and Z& re- 
ceive their major contribution from the leg tad- 
pole diagrams, in Zy 2 leg and operator tadpole 
diagrams cancel, in Zy 3 and Z& 2 the net contri- 
bution is (approximately) one power of uo, and 
in Zy 4 two powers. For Z\ yc ,Z aQ , and Zy 2 the 
agreement with improved (the dashed line) and 
even standard perturbation theory is surprisingly 
close. However, from Zy 2 onwards we find agree- 
ment only for larger momenta. It is also striking 
that the data show a somewhat steeper slope than 
the perturbative lines. (The slope originates from 
the anomalous dimensions of the operators.) This 
suggests perhaps that a more effective expansion 
parameter would be a value of a which is twice 
as large as the value used in the figures. 

For the other k values, 0.155 and 0.1515, we 
found the same results, momentum for momen- 



tum, within the error bars. One might have ex- 
pected the numbers to vary by a factor of [2,6] 
1 + (1/k — 1/k c )/2«o- This would have made a 
difference of nearly 10% between our largest and 
smallest k values. Our data do not support such 
a quark mass dependence. 

4. Conclusions 

We have seen that perturbative and non- 
perturbative renormalization constants agree 
within 20% for fi 2 > 4 and better. At the present 
value of [3 we expected the non-perturbative 
method to work for fi 2 2. 

Our data do not support taking k c = 1/8. We 
also do not see any dependence of the renormal- 
ization constants on the quark mass in the range 
of k values we have explored. 

As far as we can tell, the agreement of pertur- 
bative and non-perturbative renormalization con- 
stants is no worse than in the case of improved 
fermions [4]. Wilson fermions, however, have the 
advantage that they suffer less from statistical 
fluctuations. 
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